Mult if r act al current distribution in random diode networks 
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Recently it has been shown analytically that electric currents in a random diode network are 
distributed in a multifractal manner [0. Stenull and H. K. Janssen, Europhys. Lett. 55, 691 (2001)]. 
In the present work we investigate the multifractal properties of a random diode network at the 
critical point by numerical simulations. We analyze the currents running on a directed percolation 
cluster and confirm the field-theoretic predictions for the scaling behavior of moments of the current 
distribution. It is pointed out that a random diode network is a particularly good candidate for a 
possible experimental realization of directed percolation. 



Directed percolation (DP) is a simple model for di- 
rected connectivity in a random medium jjj, It dif- 
fers from ordinary isotropic percolation by the additional 
constraint that activity is restricted to percolate along a 
distinguished direction in space. Interpreting this direc- 
tion as a temporal degree of freedom, directed percola- 
tion may also be viewed as a reaction-diffusion process 
A <-> 2A, A — ► 0. In this process particle creation and re- 
moval compete one another, leading to a non-equilibrium 
phase transition form an active (percolating) into an ab- 
sorbing (non-percolating) phase. Because of its robust- 
ness, DP plays the role of a standard universality class 
of nonequilibrium phase transitions which may be of a 
similar importance as the Ising model in equilibrium sta- 
tistical mechanics. 

One of the simplest realizations of DP is directed bond 
percolation on a (tilted) hypercubic lattice. In this model 
neighboring sites are connected by bonds which are con- 
ducting with probability p and impermeable otherwise. 
In contrast to isotropic percolation, the bonds conduct 
only in one direction and may be realized by randomly 
distributed diodes, as shown in Fig. [I]. Obviously, the ef- 
fective conductivity of the system along the distinguished 
direction depends on the percolation probability p. If p is 
small, all cluster of connected sites are finite so that the 
conductivity vanishes in the limit of large distances. On 
the other hand, if p exceeds a certain critical threshold 
Pc, there is a finite probability to find a directed path 
between two points, leading to a finite resistance when 
averaged over many independent samples. At the criti- 
cal threshold p — p Cl the system undergoes a continuous 
phase transition where clusters of conducting paths dis- 
play a fractal structure. In this case the currents running 
on the cluster are distributed in a nontrivial way. 

The DP transition can be described by field-theoretic 
renormalization group methods ^| . In this framework 
the critical behavior near the transition can be charac- 
terized by simple scaling laws which involve three inde- 
pendent critical exponents. The order parameter in the 
percolating phase is the density of sites belonging to an 




FIG. 1: (l+l)-dimensional directed bond percolation as a 
random-diode network on a tilted square lattice. The fig- 
ure shows a realization of randomly distributed diodes with a 
conducting path from x to x' . 

infinite cluster or - using the reaction-diffusion language 
- the density of active sites. Approaching the transi- 
tion, this density scales as p ~ (p — p c ) , where j3 is the 
critical exponent associated with the order parameter. 
Similarly, the correlation length along the distinguished 
direction scales as £|| ~ \p — p c \~ v± , while the correla- 
tion length perpendicular to this direction diverges as 
~ \p — Pc\~ lJ± ■ The exponents uu and v\_ are generally 
different since introducing a distinguished direction in 
space breaks the symmetry between parallel and perpen- 
dicular degrees of freedom. In a two-loop approximation, 
these exponents are given by [|[ || 



,0 = 1-6/6- 0.01128 e 2 + 0(e 3 ) 
v\\ = 1 + e/12 + 0.02238 e 2 + 0(e 3 ) 
u± = l/2 + e/16 + 0.02110e 2 + O(e 3 ) 



(1) 



where e = 4 — d. The exponents have also been estimated 
by various numerical methods (see Table |lj) . The triplet 
of exponents (J3,i/n,v±) labels the universality class of 
<f+l dimensional DP, while other critical exponents are 
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exponent 


d = 3 


d = 2 


d = 1 


P 


0.81(1) 


0.583(4) 


0.276 486(8) 


v \\ 


1.105(5) 


1.295(6) 


1.733 847(6) 


u± 


0.581(5) 


0.733(4) 


1.096 584(4) 



TABLE I: Estimates for the critical exponents of DP. 



usually related to these standard exponents by simple 
scaling relations. This type of scaling invariance, where 
only a finite number of independent exponents is in- 
volved, is referred to as simple scaling. 

Let us now consider an electric current running on a di- 
rected percolation cluster according to Kirchhoff's laws. 
By introducing such a current the theory is extended by 
an additional physical concept with non-trivial proper- 
ties. In fact, even though the cluster itself is character- 
ized by simple scaling laws, the currents on the cluster 
turn out to be distributed in multifractal manner, mean- 
ing that all moments scale with individual exponents 
which cannot be related to the standard exponents by 
simple scaling laws. In the case of isotropic percolation, 
this phenomenon was first conjectured by Rammal et 
al. H as well as by Arcangelis et al. . In the meantime 
these results have been substantiated by field-theoretic 
renormalization group studies |9| and numerical sim- 
ulations |, 0. Recently, Stenull and Janssen fl^] 
investigated the directed case. By employing field the- 
oretic methods they determined the critical exponents 
of the moments of the current distribution. Their re- 
sults clearly indicate multifractality. The purpose of the 
present work is to verify these predictions quantitatively 
by extensive numerical simulations. To this end we ana- 
lyze a random diode network as shown in Fig. [l], which is 
a special case of the model considered in Q . Again we 
would like to stress that multifractality is not observed as 
long as one is concerned exclusively with the underlying 
DP process. An additional process, here the transport 
of electrical currents, is necessary for the emergence of 
multifractality. 

The moments of the current distribution in a d+1- 
dimensional DP process are defined as follows. Let us 
denote by x\\ the coordinate along the distinguished 
direction (=time) and by x± the vector of d perpen- 
dicular coordinates (=space). An important quantity 
in the theory of DP is the pair-connectedness function 
, x±; x'u , x'j), which is the probability to find a con- 
ducting path running from the point [x\\ , x±) to the point 
(x|| , x'j_). The union of all possible paths in a given real- 
ization of conducting and non-conducting bonds is called 
the backbone of the correlation function. Note that the 
backbone is invariant under reversal of all diodes. In the 
field-theoretic framework this invariance is known as the 
so-called rapidity reversal symmetry |l3| . 

We now assume the diodes to be ideal in the sense 
that they are perfectly isolating if the applied voltage 




FIG. 2: Example of the current distribution on a backbone 
running from x to x'. The thickness of the bonds represents 
the intensity of the current. 

is negative while they have a fixed resistance R if the 
voltage is positive. If we apply an external voltage at 
the two points x = (:E||,a;_i_) and x' = (x'^,x±) there 
will be a current running on the backbone which can be 
computed by solving Kirchhoff's laws. Obviously the cur- 
rents are distributed in a non-trivial way, especially when 
the backbone consists of many nested loops, as shown in 
Fig. ^. This current distribution has multifractal prop- 
erties which can be studied by considering the moments 

M e = Y,(h/ir, (2) 

b 

where the sum runs over all bonds b of the backbone 
while lb and / denote the current passing bond b and the 
total current, respectively. Each of these moments has a 
specific physical meaning. For example, Mq is just the 
total number of bonds of the backbone while Mi is the 
total resistance between the two points in units of R. Mi 
is the second cumulant of the resistance fluctuations and 
may be considered as a measure of the noise in a given 
realization. Finally, is the number of so-called red 
bonds which carry the full current /. These bonds play 
a special role since they are shared by all possible paths. 
At criticality the moments obey a power-law jl2| 

M^O^-xh)*/^. (3) 

Multifractal behavior manifests itself in a nonlinear de- 
pendence of the exponents ipt on their index £. The expo- 
nents ipe, just as the geometrical exponents compiled in 
Eq. (0), are known in a two- loop approximation. Thus, 
at least in principle, we can gather the e expansion for the 
exponents i^t/y\\ from the literature. This e expansion, 
however, would hardly be useful for comparison to our 
numerical data because it is valid only asymptotically in 
the vicinity of d c = 4. To produce reliable predictions for 
d = 3 or even d = 2 the e expansion has to be improved 
by incorporating additional information. Here, we carry 
out a so-called rational approximation. To be specific, 
we incorporate the rigorous feature ipi/v^d = 0) = 1 by 
supplementing the e expansion with an appropriate third 
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FIG. 3: Simulation results for the moments Me as a function 
of time in 1+1,2+1, and 3+1 dimensions. The data points 
for 2+1 and 3+1 dimensions are shifted vertically by a factor 
of 10 and 100, respectively. Each group of lines shows the 
moments Mo, . . . , M5 from top to bottom. The data points 
for M4 and M5 are almost identical. 



order term. The so obtained result reads: 



ti = 1 
v\\ 



6(4 



at 



0(e 4 



(4) 



where the ae, be and eg are ^-dependent coefficients taking 
the values listed in Table [n]. 

In order to estimate the exponents numerically, we de- 
termine the current distribution averaging over an ensem- 
ble of 5000 independent realizations of a conducting back- 
bone in 3+1, 2+1, and 1+1 dimensions. In principle the 
current distribution can be computed by solving Kirch- 
hoff's laws at all nodes of the backbone. In the case of 
isotropic percolation, where the network consists only of 
resistors, one obtains a set of linear equations which can 
be solved exactly. In the present case of a diode-resistor 
network, however, the resulting equations are nonlinear 
making it impossible to compute the currents analyti- 
cally. Thus, in order to compute the currents, we have 
to apply an iterative approximation procedure. To this 
end each site of the backbone network is connected to 
a grounded capacitor C. For a given set of voltages at 
the lattice sites, we can compute the actual difference 
between incoming and outgoing currents. This difference 
then changes the voltage of the capacitor at site i by 



dt 



hi 



* incoming -^outgoing 



(5) 



After sufficiently long time this nonlinear dynamical 
system reaches a stationary state where incoming and 
outgoing currents balance each other. Therefore, by iter- 
ating Eq. (^|), the stationary current distribution can be 
approximated accurately, mainly limited by the available 
CPU time and the precision of floating-point arithmetics. 



TABLE II: The coefficients ae, be and ce appearing in Eq. 



exponent 


d = 3 


d = 2 


d= 1 


i>o/v\\ 


1.10(2) 


1.24(3) 


1.32(3) 


i>i/v\\ 


0.95(1) 


0.88(2) 


0.82(2) 


ifa/v\\ 


0.92(1) 


0.81(3) 


0.67(3) 


i/>a/v\\ 


0.91(1) 


0.79(3) 


0.63(3) 




0.91(1) 


0.78(3) 


0.62(3) 



TABLE III: Numerical estimates of the critical exponents. 



Moreover, the convergence during the iteration proce- 
dure can be controlled easily by monitoring the quantity 
J2 b |AJj,|, which should tend to zero as the stationary 
state is approached. Using Newton's iteration scheme, 
this quantity actually decreases until it reaches a value 
of about 10 -10 , which is essentially determined by the 
limited machine precision. At this point we stop the it- 
eration procedure and compute the moments Me, assum- 
ing that the errors caused by the iteration procedure are 
much smaller than the statistical and finite-size errors 
generated by the DP process itself. 

In the simulations we observe that the average size of 
the backbone increases rapidly with x'^ — x\\, especially 
in 3+1 dimensions. While the numerical effort for gen- 
erating conducting backbones between two given points 
plays a minor role, the algorithm is mainly limited by 
the CPU time needed for the dynamical approximation 
and by the memory capacity for storing the configura- 
tion of the backbone, the currents, and the voltages. Be- 
cause of these limitations, we could only go up to dis- 
tances of x'n — a; || < 75 in 3+1, 2+1 and 1+1 dimen- 
sions. Clearly, in this limited range corrections to scaling 
still play an important role. Nevertheless we find that 
the moments exhibit very clean power laws (see Fig. |^), 
which allows us to estimate the critical exponents ipg i n 



the limit x'„ 



'II 



00 by standard extrapolation tech- 



niques. These estimates are summarized in Table III. 

Figure [| shows the field-theoretic estimates as given in 
Eq. (||) compared to the numerical results. The error bars 
include the statistical error of the simulation combined 
with the expected error of the extrapolation procedure. 
In 3+1 dimensions we find excellent agreement. The nu- 
merical values and the field theoretic predictions almost 
coincide. As expected, the agreement is less pronounced 
in 2+1 dimensions. The analytic estimates lie slightly 
outside the error bars of the simulations. Though not 
quantitatively correct, field theory yet predicts the right 
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FIG. 4: Comparison of simulation and field theory for d = 2 
(top) and d = 3 (bottom). The field theoretic results are 
indicated by squares. The numerical ones are symbolized by 
diamonds with error bars. 

shape of the dependence of ipi/v\\ orL ^- I n 1+1 dimen- 
sion, the field theory fails to predict the numerical values 
of the multifractal exponents accurately. We therefore 
refrain from showing the respective graph. All in all, our 
simulations support the validity of the field theoretic re- 
sults, at least for dimensions to that the e expansion can 
be continued reliably. Our simulations clearly confirm 
the existence of multifractality in current distributions 
running on DP clusters. 

As in the field-theoretic case, we can check the consis- 
tency of our results. For example, as shown in fulf , the 
zeroth moment in a (d+l)-dimensional system, which is 
just the mass of the backbone, can be expressed in terms 
of the backbone dimension Db = d + 1 — 2/3/i/± (cf. 
Ref. @) by 



V>o = v±(D B -l + z). 



(0) 



where z = un/v±. Using the values listed in Table |, we 
obtain for ipo/^\\ the values 1.11 in three, 1.23 in two, and 
1.307 in one dimension, respectively. These values are in 
perfect agreement with the estimates listed in Table HI, 

To summarize, we have studied the moments of cur- 
rents running in a random diode network at the critical 
point by numerical simulations. Our results fully confirm 
recent field-theoretic results, supporting the claim that 
currents running on a critical DP cluster indeed show 
multifractal behavior. 

We would like to point out that a random diode net- 
work might be a good candidate for an experimental re- 
alization of directed percolation. In particular the total 
resistance Mi is easy to measure. In an experimental 



realization various idealized assumptions of the model 
shown in Fig. |l| can be relaxed. For example, not all 
diodes need to be polarized along the distinguished di- 
rection, rather it is sufficient for them to be preferen- 
tially oriented in a statistical sense. In addition, the re- 
sistance of conducting diodes may be stochastically dis- 
tributed. Nevertheless the critical exponents of such a 
system at the transition are expected to be the same as 
in DP. Even inhomogeneities, which cannot be avoided 
in any experiment, are not expected to play a role since 
they would result into an annealed noise which is irrel- 
evant under renormalization group transformations (see 
Ref. fl5[] ). Because of these favorable properties, random 
diode networks are particularly suitable for a possible ex- 
perimental realization of directed percolation. It should 
be emphasized that so far no experiment is known where 
the critical exponents of directed percolation could be 
verified in a reliable way Jjl], |l6). Therefore, the experi- 
mental realization of a random diode network would be 
a particularly challenging task. 
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